A COMMON FORMULA FOR CERTAIN 
GENERALIZED HANKEL TRANSFORMS 



MARIO GARCIA ARMAS 



Abstract. We study the generalized Hankel transform of the family of se- 
quences satisfying the recurrence relation a„+i = (a + „^^ )cin,- We apply 
the obtained formula to several particular important sequences. Incidentally, 
we find a connection between some well known formulas that had previously 
arisen in literature in dissimilar settings. Additionally, given a non-zero se- 
quence satisfying the above recurrence, we evaluate the Hankel transform of 
the sequence of its reciprocals. 



1. Introduction 

We recall some terminology from the theory of Hankel matrices. Given a se- 
quence (a„)^o, we consider the doubly-indexed sequence of Hankel matrices Hn^\ 
n = 1, 2, . . .; A; = 0, 1, . . ., defined by 



afc+2 



afc+2 



flfe+a 
afc+4 



Ofe+ri 
Ifc+n+l 



(1.1) 



\ Ofc+Ti-l Ofc+n ak+n+1 ■ ■ ■ afc+2n-2 / 

The (generalized) Hankel transform of the sequence (an)5JLo the doubly-indexed 



sequence of determinants dn ' = detiJn (for a similar treatment, see [T^lUnj)- It 
is important to mention that several authors refer to the Hankel transform only as 
the sequence dn = detHn^ (see, for example, [3l [T7]). 

The theory of Hankel matrices have beautiful connections with many areas of 
mathematics, physics and computer science (see, for example, [T^ HOI [H] ) . Al- 
though Hankel determinants associated to sequences had been previously studied 
(see, for example, [T]), the term Hankel transform was introduced and first studied 
by J. W. Layman in [T7]. Several later studies of Hankel transforms of sequences 
have appeared in literature (see, for example, [3l [71 [9l [TOl [TTl [T9] ) . 

In the evaluation of Hankel determinants, several techniques have proved to be 
useful. For an extended set of tools, as well as a significant bibliography, please 
refer to Sec. 2.7], [HI Sec. 5.4] and [H]. 

In this paper, we study the generalized Hankel transform of a sequence (a„)^Q 
satisfying 



(3 



■7 



Vn > 0, 



(1.2) 
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for some complex nmiibers a, (3, 7. 

Moreover, we study particular choices of these values, for which important se- 
quences are obtained. Some of these Hankel transforms have been previously stud- 
ied in quite independent settings. Consequently, one of our main goals is to illustrate 
the hidden connection between them. 

In addition, we note that the sequence of reciprocals of any non-zero sequence 
satisfying (|1.2p is also defined by a similar recurrence relation for some parameters 
depending on a, /3, 7 (see Equation (|3.27p ) and consequently, we evaluate the Hankel 
transform of the sequence of reciprocals. 

It is worth mentioning that in general, it cannot be obtained a closed product 
form evaluation for Hankel transforms of sequences satisfying the recurrence 

q{n) 

where p and q are arbitrary polynomials with maxjdegp, deg g} > 1. The non- 
existence of product evaluations can be proved by means of the procedure described 
in [TU|, i.e., computing the first terms of the transforms associated to integer (or 
rational) sequences and factor them into primes. The presence of extremely large 
primes in the factorizations shows that there cannot be a product form evaluation 
of the transforms. For example, the Hankel transform of the sequence a„ = ('^"^^) 
does not have a product form evaluation (see |10j ) and however, it is defined by the 
recurrence 

3(3n + 2)(3n + 4) ,^ ^, 

2(n -I- l)(2n -|- 3) 

Analogously, it can be proved for example, that the Hankel transforms of the se- 
quences defined by the recurrence relations 

a„+i = — a„ (1.5) 



and 



1 



(1.6) 



cannot be evaluated in closed product form. 

Therefore, it is unlikely that general closed products form evaluations of Hankel 
transforms of sequences satisfying (|1.3p could be obtained outside from the linear 
case (i.e., degp = degg = 1), which is the one studied in our paper. 

2. Basic Computations 

Throughout this paper, we consider a finite product Jls=a Cs = 1 when b < a. 

Let (a„)^Q be a sequence satisfying (|1.2I) . For simplicity, we may assume that 
uq — 1. It is quite obvious that the results we derive in this section can be straight- 
forwardly extended to the general case (see Remark 12.41 below). 

Note that for every to > n > we can write 

" f ais + ^-l)+f3 \ "fr7 a(. + n + 7-l) + /3 \ 
a,„ = a„ M =«" M — — ■ (2.1) 

Let us consider the Hankel matrix i?!*^'' — {'mi.j)i<,ij<,n, with the notation of 
p.ip . Note that we have to^ — ai+j+k-2 and consequently, from (|2.ip we obtain 
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the relation 

(k) 

We now obtain a first tentative expression for c?„ . 
Lemma 2.1. 

n-1 

[a(*-l)-/3]""' 

d-n^ = afeOfc+i • ■ • afc+n-i — • (2.3) 

Y[[i + k + -f-lY [2n-i + k + 'y~2y 
1=1 

Proof. This formula could be derived, after some work, from the generalization of 
the Cauchy determinant described in [6, Thm. 4.3]. 

Alternatively, a different (and more elementary) approach will be outlined. Let 
M„(a, b, c) denote the n x n matrix with entries 

< = n "^'^i^\ ^,J = l,2,...,n; (2.4) 

s + J + c 

where a, b, c are complex numbers. Then, it can be proved that the following 
recurrence relation holds: 

(n - W (nr - h)^^^ 

detM„(a,6,c) = ^^-^ '-^ det Af„_i(a, a + &, c + 2). (2.5) 

Y[{i + l + c){i + 2 + c) 
1=1 

This can be achieved after the following sequence of steps: 

• Subtracting the (fc — l)-st column of A/„(a, b, c) from the fc-th one for k — 
n,n — 1, . . . ,2, in that order. 

• Expanding by the minors of the first row. 

• Taking common factors from rows and columns out of the determinant of 
the lower right (n — 1) x (n — 1) block (These common factors form the 
fraction on the left of the RHS of Equation (|2.5p ). 

• Noticing that the resulting matrix is M„_i(a, a + b,c + 2). 

As a consequence of Equation (12. 5p , we obtain the formula 

n-1 

[a{i-l)+ac-b]"~" 

detM.n{a,b,c) = ^^ . (2.6) 

Jl + 1 + c] " [2n - i + c] ' 



The proposition follows easily from Equation (|2.6p after the substitution 

(a, b, c) := {a,aj + ak~2a + /3,-f + k~ 2). (2.7) 

□ 

(k) 

Let us find a nicer expression for dn . We make the following observation: 



4 



MARIO GARCIA ARMAS 



Remark 2.2. It is easily seen that dn = implies c?„ =0. In fact, the numerator 



of the fraction on the right of (|2.3p does not depend on k and on the other hand, 
if Gi — for some i G {k, . . . ,k + n — 1}, then Equation (|1.2p yields aj = for aU 
3 > i- 

Suppose now that we have S/^ ^ for some j > 0. Using Equations p.ip and 
together, we compute the ratio 



0"j-]-n 2—1 
— X 



n + 7- 1]' [2" -« + J +7-2]' 



— X 



]J[i + j+7]' [2r^-^+J+7- 1]' 

i=l 

n-1 



J 

1=1 

n n— 1 

[][a(z + j+7-l)+/3] llii + J+l-l) 



i=l 



1=1 

n Q;(i+j + 7- l) + /3 
i + j + 7 + n-2 ' 



Joining the last result with Remark 12.21 we conclude that the relation 

= n (2.9) 

is always valid and consequently, we obtain the formula 

7=0 r= 1 i + + 
Now, let us focus on dfi ^ . Following Equation (j2.3[) . we can write 

n-l 

Kz-l)-/3]"-^ 

4°) = anai • • • ^n-i ^^ZT^ ■ (2-11) 

[|[i + 7-l]' [2n-z + 7-2]' 

z=l 
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It is easy to see that, similarly to Remark [221 — implies d^n+i — ^- Again, 
if we have cff^ / for some j > 1, we compute the ratio 



ajj\ X 



i=l 

\[[a{i-l)-p\ 

i=l 



f\\a{^-l)-^3\ 



i=l 



., A a(z + 7 - 1) + P ^ ^=l 

i + 7- 1 J 

1=1 

y [i + 7-l][«+.?+7-2][i + j+7-l]' 



Hence, we can conclude that the relation 



[* + 7-l][* + J+7-2][z + j+7-l] ^ ^ 



is always valid and therefore, from the obvious d!f^ = oq = 1 we obtain 

^(0) ^ TT i[a(i + 7-l)+/3] 

" i<.ii-i [* + 7-l][*+J+7-2][z+j+7-l]' 
We summarize the obtained results in the following theorem. 



(2.14) 



Theorem 2.3. Let (a„)5^Q be a sequence such that — 1 and satisfies Equation 
(|1.2[) . Then, its generalized Hankel transform is given by 

^(k) ^ -Q z[a(z + 7-l)+/3][a(*-l)-/3] 



_i [* + 7-l] [z+j + 7-2] [z+i + 7-l] 



l<i<j<n 

/y=j^A a(z+,+7-l) + /? _ 

i+7+7 + n-2 ^ ^ 

Remark 2.4. If we drop the restriction oq — 1, then it is easily seen that the 
expression for d^"^ will be multiplied by Oq . 

3. Applications to Particular Sequences 

We now consider several applications of Theorem 12.31 to particular important 
sequences. In some cases, we show how Formula (|2.15p applies to derive some 
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well known closed product form evaluations, which have been obtained before us- 
ing several independent methods. We emphasize that empty products are always 
considered to be 1. 

Application 1. We start with the sequence a„ = l/{n + 1), whose associated 
matrices -ffn'^'' are (generalized) Hilbert matrices. This sequence satisfies (|1.2p for 
a=l,/3=— 1,7 = 2 and thus, its Hankel transform is given by 



.o k—1 n 

1 TT TT * + 



4'^= TT T — T^Ar rxTTTT- ■ 



After some elementary transformations, we obtain the identities 

ra-l 

n 



l<j<j<«-l L ' J L ' -'J L ' J ' J TT(*+i^l) 



k—l n . , . n . . 

* + J TT * + .7 ^ 1 



TTTT ^ = TT 



which allow us to write d„ in the more familiar form 



n-l 

\2 



(3.1) 



= ^^^^^ ■ (3.2) 

n (i+j+fc-i) 



This is a very well known formula and can be proved by several methods (for 
some historical remarks, see 18, vol. Ill, pp. 311]). For example, it can be easily 
derived from Cauchy's double alternant (see |15l Eq. 2.7]). An extensive literature 
exists on Hilbert matrices and their generalizations. For an interesting compilation 
of results about Hilbert matrices, please refer to [4j. 

More generally, it is easy to see that the Hankel transform of the sequence 
a-n = l/(An + fi) is given by 



n 



l<'i<j <n 



+ k] 

-1 [« + k] +j + - 1] [i+j + n] 



k — l n 

nn 

3=0 »=1 



i + j + K - I 
z + + l' 



(3.3) 
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where n = n/X (see Remark [2^ . Again, using the identities 

11(^0 

l<i<j< 



n-1 

2 



n-i[^ + -][^+^+--^][^+J + -] ^(^+, + .~2)' 

fc-1 ri . . , n „ 

n-p-r I + J + K - I _ -p-r t + J + K — 2 

. i + i + n + K.-l ^ i-l- i + J +^ + ^-2 ' 

J=0 1=1 «J = 1 

we obtain the more general (and well known) formula 

n-1 

4!"^ - ^ X , (3.4) 

which could also be derived, for example, by means of the partial fractions method 
(cf. Thm. 4.1]). 

Application 2. Let a„ — 2/(n^ + 3n + 2) be the sequence of the reciprocals 
of triangular numbers shifted by one. It satisfies (|1.2p for a — 1, P — ~2, j — 3. 
Thus, its Hankel transform is given by 



■2 r ■ i 1 1 k — l n . , . 

rf(fc) ^ n I + X TT TT ' + ^ 

1<1<J<TI-1 L J L J L J j=0 4=1 

By considering the identities 

n-1 



(3.5) 



TT 6 ' [i + 1] 

J- i r,- j_ 9i r,- ^ o ^ 1 1 li J- 1 J- 9i ~ 



i=l 



1^2 



i<i-i-i[^ + 2][. + , + l][^+,+2] j!^^^^^.^ 



t + J f n + k\ -f-r i + J 



j=0 i=l ' ' ' \ ^ ij = l •' 



n n 1 4^ n + 1 ( n ) n 

3=0 i=l 

we obtain the simpler formula 



n-1 

^2 



'^"'='"rn ) "^^^^^ — • 

We could have easily evaluated the Hankel transform of the slightly more general 
sequence a„ = (n + A)^^(n + A + 1)^^ by the same method. However, note that 
the quadratic polynomials in the denominators are carefully chosen such that the 
sequence satisfy (jl.2p for some a, /?, 7. In fact, if we consider the apparently similar 
sequence a„ = (n+l)^^^ -^q note that its Hankel transform cannot even be evaluated 
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in closed product form! This can be proven by the techniques from [10] (see also 
Introduction) by noting for example the factorizations 

^(0) ^ 647 

28 ■ 36 • 52 



^3 - 98 . C.6 . t;2 ' (^■'^) 



(0) _ 179-179357 

"5 - 220.36.510.75' 

(0) _ 23 . 1280587616051046200369 
^ - 236 . 322 . 510 . 714 . . 132 ■ ('^■9) 

The amazingly large primes in the numerators confirm our claim. 

Application 3. Consider the sequence a„ = 1/nl. It satisfies (11.21) for a — 0, 

/3 = 1, 7 = 1 and therefore, we have 

fc-i 



4'=) = (-i)(s) rr — - X TT fr — ^ — - . 



(3.10) 



l<i<j<n-l 

Alternatively, the following identities hold: 



1 i! 

n u+j] " n ' (3.11) 

l<i<j<n-l I ' -J J L ' Jj i=0 ^ ' 

t^" 1 _ (z + n-1)! 

-^i-^^z + j+Ji-l ~ lMz + fc + n-1)!' ^^'"^^^ 

Hence, we obtain 

n — 1 .| 

= (-i)(;) n r-ri — ' (3.13) 

i—O 

which recovers the formula from [2J Thm. 1.3]. 

Application 4. Let a„ = (n + 1)~^(^) be the sequence of Catalan numbers. 
It satisfies (|1.2p for a — 4:, P — —6 and 7 = 2. Therefore, its Hankel transform is 
given by 

Now, the left product reduces to 1 and the right product can be rewritten as 

n ("5) 

i<«<i<fc-i 

for /e > 0, which is obvious from the identity 

i=l 4=1 

Accordingly, we obtain the well known formula 

4'=^= n '-^^^ (3-17) 

l<i<j<k-l •' 

which was primarily found by Desainte-Catherine and Viennot ([8]), who also gave 
a combinatorial interpretation for this transform (see also |14j for further general- 
izations). It is also proved in [20], by means of the Dodgson's condensation method 
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( [151 Prop. 10]) and discussed in [16, Thm. 33], with some additional remarks. The 
cases k = 0,1 and 2 are also studied by Aigner in [1], where the author describes 
a beautiful generalization of Catalan numbers (called Catalan - like numbers) in- 
spired by the property dn'' = dn"^ = 1. 

Application 5. Let a„ = {^^^ be the sequence of even central binomial coeffi- 
cients. It is immediate to see that (a„) satisfies (|1.2|) for a = 4, /3 = —2 and 7=1. 
Hence, its Hankel transform is given by 

1<J<J<»1— 1 L J L -J 



4(^ + j)-2 
i + j + n — 1 



(3.18) 



The left product is readily seen to be 2" ^ 
rewritten as 

i + 3-1 



3=0 1=1 

As for the right product, it can be 



2x n 

l<i<j<k-l 



2n 



for fc > 1. This can be deduced directly from the identities 



n 



4(z + i)-2 
i + j + n — I 



n 



i+j-l + 2n 



n 



4i - 2 
i + n — I 



= 2. 



Thus, we are able to obtain the formula 

2"~\ 



i+j-l + 2n 



i + j-1 



if = 0; 
if fc > 1. 



(3.19) 

(3.20) 
(3.21) 

(3.22) 



2" X Y[ 

l<i<j<k-l 

2m\ _ nf2m-i\ j-gcover the formula for the Hankel 



Taking into account that (^) = 2( 
transform of odd central binomial coefficients from [20l Eq. 1.5] (see also [121 Eq. 
13]). 

It is worth mentioning that, given the identity (^^) = (—1)" 2^" ( l^^), the 
Hankel transform could also be computed using 15, Eq. 3.12]. For interesting 
connections of this Hankel transform with combinatorics and algebra, see [T] and 
[l2] respectively. 

Application 6. Let a„ = (^), where A G C. It satisfies (|1.2p for a = — 1, 

(3 = X + 1 and 7 = 1; consequently, its Hankel transform is given by 



l<i<j<n-l L' 



-A-1 



A] 



nn 



^ + J - A - 1 
i + j + n - 1 



(3.23) 



This evaluation, as well as the one from the next example, could be also derived 
after some work from [15[ Thm. 26], where the author describes a powerful set of 
tools for evaluating determinants of matrices with entries given by binomial coeffi- 
cients. Some additional results in binomial determinant evaluations can be found 
in [51 [6]. It is also worth mentioning 13], where the authors describe interesting 
combinatorial interpretations of a wide class of binomial determinants. 

Application 7. Let a„ = ("m^)' "^^^re to is a nonnegative integer and A is 
either an integer with A > m or A G C\Z. Note that oq — (^) and (a„) satisfies 
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the relation 

fln+i = ( 1 + , , ^) Vn>0. (3.24) 

Thus, by applying Theorem 12 . 31 together with Remark l2.4l we easily see that the 
Hankel transform of (a„) is given by 



dW^f^Y ]J i[i + X][i-l-m] 



V- -/ i<.<j<«-i [* + A - m] [z + j + A - m - 1] + J + A - m] 

X TTrr i±i±^ . (3.25) 

j=O i=l ' ' 

Then, it is straightforward to note that di^^ = if (and only if) n > to + 2. 
Another interesting consequence is that in the limiting case n = m + 1, we have 

rfL'li = (-l)("^'^ (3.26) 

which can be proved without difficulty after some manipulations (note that the 
second product reduces to 1). 



3.1. Hankel Transforms of Reciprocal Sequences. Wc finish this section by 
noting the following beautiful property: if (a„) is a non-zero sequence satisfying 
Equation (|1.2p with a ^ 0, then it is obvious that the reciprocal sequence (a^^) 
satisfies the relation 

- + 7 Wl ^ \ 1 ^3.27) 



a„+i an + + (3 an \a n + ^f + ^^J a. 

Corollary 3.1. Let (a„)5^o non-zero sequence such that qq — I and satisfies 

Equation (|1.2p for some a =/= 0. Then, the generalized Hankel transform dii^ of 
("n^)r=o diven by 

-j-j- i[i + -f-l][a{i-l) + (3] 



l<^<J<n-l + 7 - 1) + /3] [ait +J+J-2) + 13] [a{i + j + 7 - 1) + /?] 

k — 1 n 

V TTTT 

a{i + j + 7 + n - 2) + /3 



k-l n . . 

xTTTT . '^'o^ (3-28) 

11 11 ^("i ^ T _|- ^ ^ T7 - 9"! ^ « ^ ^ 

3=0 i=l 

Proof. The result follows from Equation (|3.27p and Theorem 12.31 □ 

Remark 3.2. The formula remains valid in the case a = 0; an easy way to see this 
is by making a — s- in Equation (j3.28p . 



As an immediate consequence of Corollary 13. 1[ we evaluate some non-trivial 
generalized Hankel transforms: 



A COMMON FORMULA FOR HANKEL TRANSFORMS 



11 



• Let Gn = {n+ 1) (^^) ^ be the sequence of the reciprocals of Catalan num- 
bers. Then, its Hankel transform is given by 

M) _ 1 n ^[^ + l][2^-5] 

" On(n+fe-l) ii 



[2i - 3 


] 


[2i-l] 


[2(i + i)-3] 


[2(i + i)-i] 



2„(„+fe-i) . [2i - 1] [2(i + i) - 3] [2{i + j) - 1] 

fc-l n -I'll 

xTTTT o.- ! o - (3-29) 

• Let a„ = (^!^) ^ be the sequence of the reciprocals of even central binomial 
coefBcients. Then, its Hankel transform is given by 

rf(fe) = I TT 

" 2n(n+k-l) 11 

l<z<j<n- 

• Let a-a = (^) , with A G C\N. Then, its Hankel transform is given by 
"^^"^ iJ-Li [^-A-l][z+,-A-2][^ + ,-A-l] 

fe— 1 n . . 

xTTTT ^ . (3.31) 

i + j + n- X-2 ^ ' 

Acknowledgments. The author would like to thank Prof. B. A. Sethuraman for 
introducing him to the subject and for so many fruitful comments and discussions. 
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